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Abstract 

Topologically stable cellular partitions of D dimensional spaces are studied. A complete 
statistical description of the average structural properties of such partition is given in term 
of a sequence o/ -j — 1 (or ^^j^) variables for D even (or odd). These variables are the 
average coordination numbers of the 2k-dimensional polytopes (2k < D) which make the 
cellular structure. A procedure to built D dimensional space partitions trough cell-division 
and cell-coalescence transformations is presented. Classes of structures which are invariant 
under these transformations are found and the average properties of such structures are 
illustrated. Homogeneous partitions are constructed and compared with the known structures 
obtained by Voronoi' partitions and sphere packings in high dimensions. 

1 Introduction 

We study topologically-stable division of any dimensional space by cells. Such systems have 
minimal incidence numbers. Configurations with higher incidence numbers are topologically 
unstable because they can be splited into configurations with the minimal incidence numbers 
by infinitesimal local transformations. In the literature these cellular partitions are known 
as "froths" since in two and three dimensions {2D and 3D) the soap froth is the archetype 
of such structures. A 2D froth is a space-filling cellular partition made of irregular polygons 
where on each vertex are incident three polygons. A 3D froth is a polyhedral partition of 
space where on each vertex are incident 4 polyhedra. In general, a D-dimensional froth is a 
partition of space in irregular polytopes, where on each vertex D polytopes are incident. 
Cellular structures with minimal incidence numbers always appear when the space is filled 
by cells without following any special symmetry. Therefore, froths are the typical structures 
of any disordered partition of the space in cells. 

A broad class of disordered natural and artificial cellular systems have the topological 
structure of froths |]T], |^, ^. Examples in two dimensions are magnetic domains in garnets 
films, Bernard-Marangoni cells in thermal convections, biological tissues, cuts of polycrys- 
talline metals and ceramics, emulsions, the subdivision of territory in administrative regions 
or in national states, geological structures and 2D soap froth (which is obtained by squeezing 
the soap foam between two plates) 0, il |]. In three dimensions, examples are biological 
cells, poly crystalline metals and ceramics, foams @, 0, H. Moreover, the structure of any 
packing (of hard spheres or atoms, for example) is the dual of a cellular system (which 
can be generated, for instance, by using the Voronoi construction ^ around the centers of 
the packed elements). In general, cellular systems generated by packings elements without 
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the use of any specific symmetry have structures which are topologically froths. It follows 
that, among the examples of 3D froths one can include amorphous metals, glasses and some 
crystalline structures such as the tetrahedrally closed-packed phases (t.c.p.) [|T0|, |11|. 

Froths in spaces with dimensionality higher than 3 are relevant in information theory 
and signal processing |T^, [TB|. Indeed, an information can be associated with a point in 
an A^-dimensional space. To transmit and recover the information in presence of noise one 
must put the points in the A^-dimensional space, separated by a certain distance which must 
be larger than the additional noise. Therefore to each point (information) is associated a 
finite volume and the entire space is subdivided in cells each one containing one encoded 
information [|l^. The energy necessary to transmit an information is proportional to the 



distance of the representing point respect to the origin. An efficient coding, which minimizes 
the energy, organizes the volumes associated with the different information in the closest 



possible packing of similar cells around the origin |T3 



Dense packings of equal cells in high dimensions have also applications in the study of 
analogue-digital converters. In this case, the space of the continuous analogical variables is 
quantized in a system of cells and the volume inside each cell is associated with one digital 
information. The quantization error is associated with the extension of the interface between 



the cells and with the distance between the center of a cell and its vertices 



High dimensional partition of space has also application in neural networks and complex 



system dynamics fl^, |15|, |T^, |T8|. Some relevant properties (such as the storage capacity 
in neural network and the slow aging dynamics in glasses) are associated with the subdivision 
of the phase-space in basin of attraction around the stored information or the minima of the 
energy. 

Two and four dimensional froths (and their duals: triangulations and simplicial decompo- 
sitions) have relevance in quantum gravity |T^, ^ |2^, Q . Here the continuous space 



is divided into cells and the functional integration over all equivalence classes of metrics is 
replaced with a summation over all the triangulations of the given manifold. 

Despite the broad variety of system which are topologically froths and the large amount 
of studies devoted to them in the hterature, very little is known about the structure of froths 
in dimensions larger than D = ?>. Froths are disordered cellular structures where the cells 
are highly correlated. These correlations essentially come from the space-filling condition 
which locally constraints the cells to pack without leaving any empty space, and globally 
constraints the froth to tile a manifold with a given curvature. In this paper we study how 
these local and global conditions determinate the average topological properties of the froth 
structure and we construct froths in any dimension by cell-division and coalescence trans- 
formations. The aim of the present paper is to investigate the average structural properties 
of classes of homogeneous partitions of high dimensional Euclidean spaces and to give ana- 
lytical instruments and methodologys for the the investigation of the topological structure 
of froths in spaces of arbitrary dimensions and curvature. 

The plan of the paper is the following: In Section |^, the hierarchical organization of 
topologically stable divisions of space in cells is studied; In Section ^, we discuss a way to 
generate or modify D-dimensional froths by cell-division and cell-coalescence transforma- 
tions; In Section ^ the fixed points of such transformations are studied and the properties 
of the associated structures are illustrated; The construction of homogeneous partitions and 
the comparison of their properties with known structures, is done in Section |^. 
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2 Hierarchy in the cellular structure 



A froth in an arbitrary dimension D is a cellular structure where the incidence numbers (i.e. 
the average number of elements which are incident on a given lower-dimensional element) 
are fixed by the stability condition at the minimal value. The cells of a froth in dimension D 
are D-dimensional irregular polytopes packed together to fill space. The boundaries of these 
cells are made with (D — l)-dimensional polytopes which are bounded by (Z^ — 2)-dimensional 
polytopes and so on up to the zero-dimensional elements which are the vertices. For example, 
a three dimensional froths is made with 3D polyhedra (the cells) which are bounded by 2D 
polygons (the faces) which are bounded by ID elements (the edges) finally bounded by QD 
elements (the vertices). A characterization of this D-dimensional structure can be given in 
term of the numbers of D-dimensional polytopes which are making the froth, in term of the 
average numbers of {D — 1) -dimensional polytopes making the boundary of a given cell and 
so on counting the number of polytopes making the boundary of the boundary, etc. 

The boundary of any /c-dimensional polytope of the froth is also a froth in a (A; — 1)- 
dimensional elliptic space. The D-dimensional froth is therefore a graded topological set: it 
contains £)-polytopes, the cells, which are tiling a space which can be Euclidean, elliptic or 
hyperbolic. The boundary of each cell is an elliptic D — 1 surface which is tiled by a {D — 1)- 
dimensional froth, whose cells are {D — l)-polytopes which are the interfaces bounding the 
original cell and separating it from its topological neighbours. Each interface of the [D — 1)- 
froth is an elliptic {D — 2)-froth of {D — 2)-polytopes, which are separating the cells from 
their neighbours. The graded topological set terminates with edges (segments or convex 
1-polytopes), bounded by 2 vertices (or 0-polytopes). 

Let denote with Ck the number of A;-dimcnsional cells in the froth ( Co number of vertices, 
Ci number of edges, C2 number of faces, C3 number of polyhedra ... Co number of D- 
dimensional cells.). Let denote with (rijj) (for i < j) the average number of i-dimensional 
cells which arc surrounding and making the boundary of a j dimensional cell ((^0,1) — 2 
number of vertices surrounding an edge, (^1^2) number of edges per faces, etc.). 

The average froth structure is characterized by the numbers of polytopes Cj and by the 
valences (riij) (with < i < j < D), which are therefore the variables of the problem. 
The total number of these variables is ^{D + 1){D + 2), but they are related by the Euler 
equations and constrained by the stability condition. In particular, the numbers of elements 
in the froth Ci are related through the Euler relation: 



where xd is the Euler-Poincare characteristic which is associated with the space curvature. 
For D even, opposite signs of xd correspond to spaces with opposite curvature: xd > 
corresponds to an elliptic D-dimensional space and xd < corresponds to an hyperbolic 
space. On the other hand, for D odd, this relation between the sign of xd and the space 
curvature does not holds anymore. 

An Euler relation is also satisfied for each froth of the graded topological set. That gives 
a set of relations for the quantities {riij) 



D 



(1) 



i=0 



J-1 



y: i-iTM = 1 - (-1)^ (with J = 1, 2, D). 



(2) 



i=0 
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Here the factor 1 — (— l)"' = ■ySj^^^^^'^^ is the Euler-Poincare characteristic for the surface of 
a J-dimensional sphere (which is a (J — l)-dimensional elhptic space). 

The numbers Cj and the averages {uij) are related by the stabihty condition 

{n,^j)Cj (with i < J < (3) 

where rii^i = 1. The binomial coefficient in the left-hand side of Eq.(^ is an incidence number 
(number of j-dimensional polytopes incident on an z-dimensional polytope) which is fixed 
at the minimum value by the stability condition (there are {D + 1) edges and ('^2^^) faces 
incident on each vertex, D faces incident on each edge, etc.). 

By contrast, the coordination numbers (^i.j) with i < j are variables (except for no.i: 
every edge is bounded by two vertices and therefore no,i = 2). These variables are not all 
independent and their range of variability is severely restricted by relations (|I]), (0) and (^) 
(for example, in a 2D infinite Euclidean froth we have (^1,2) = (^^0.2) = 6 as consequence of 
the Euler relation). 

One can show that p5| a complete topological characterization of the average structure of 



a D-froth is given by a set of y — 1 (or ^j^) for D even (or odd), independent variables: the 
even "valences" X21 = {n2i-i,2i) with 21 < D. These valences are the average coordination 
numbers (average number of neighbours) of the 2/-dimensional polytopes in the froth. These 
are free variables. On the contrary, the coordination numbers for the odd dimensional 
polytopes (the odd valences) are given in terms of the even valences by the relations 

^i^2jyr^-^5^ + ----(-l)-'^^ + (-l)% = l-(-l/ (for./=l,2,....D). 

(4) 

(Which is obtained from Eq.(^ associated with Eq.(^ and by using the definition = 
{nk-i,k)-) For J = 21 + 1 odd, the left hand term in Eq.(^ is equal to 2 and Eq.(^) fixes the 
value of the odd valence X21+1 in terms of the even ones. When J = 21 even, the left hand 
term is zero and therefore the even valences X21 are free variables. 

In k dimensions the polytope with minimal coordination number is a simplex with k + 1 
neighbours. Therefore, the valences must stay into the range + 1 < X^ < 00 ■ The 
average structure of a D-dimensional froth is characterized by a sequence of even free va- 
lences {X2, X4, Xq }. For any given sequence of even valences {X21} the odd ones can 

be calculated by using equations (^). But only sequences which generate odd valences with 
X21+1 >2l + 2 are admissible. This condition strongly constraint the accessible values of the 
even valences. 

The Euler's formula (|ip, associated with Eq.(^ and Eq.(^), gives an additional relation 
between valences 

In even dimensional spaces, the sign of the Euler-Poincare characteristic xd is associated 
with the space-curvature. The sign of the term inside the brackets in the left-hand side of 
Eq.(^ is the same of xd (because Cd > 0). Therefore, any two regions of the valences' space 
which have different signs of the bracket term (i.e. of the xd) correspond to two froths on 
manifolds of opposite Gaussian curvature. For example, in the 2D case (where the average 
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structure of the froths is described by the parameter X2) equation (|^) gives 



(6-X2)^ = X2 , (6) 

which indicates that 2D froths with X2 < 6 are tihng elhptic surfaces, whereas froths with 
X2 > 6 are tihng hyperbohc surfaces (see fig. (la)). In four dimensions Eq.(^ gives 

(where we used Eq.(H) to express X3 as a function of X2). Equation (|^) indicates that 
the region in the parameter space {X2,X4} below the hne X4 = 5|5^ is associated with 
AD froths which are tihng elhptic manifolds (x4 > 0), whereas the region above this line 
correspond to froths tiling hyperbolic manifolds (x4 < 0) (see fig.(lb)). 

Note that Eq.(^ is a constraint on the sequence {X21} given by local conditions (it 
concerns the average properties of a 2/ + 1 dimensional polytope in terms of the properties 
of the lower dimensional elements that are making it). Whereas, Eq.(^ is a constraint on 
{X21} given by the global curvature of the manifold that the froth is tiling. 

Note also that relation (|^ has a singularity in X2 = 5. This is associated with the 
existence of polytope {5,3,..} (Schlalfly symbols [^) up to = 4 only, as explained in 
Appendix 0. 



3 Cell-division and cell-coalescence transformations 

In this section we build D- dimensional froths by using the cell-division transformation and 
its inverse (the cell-coalescence). This is a local transformation that changes the structure 
of the froth but leave unchanged the global topological properties (the curvature of the 
manifold tiled by the froth or -equivalently- the parameter xd)- By using cell division and 
coalescence it is therefore possible to generate different froths which are tiling topologically 
identical manifolds. 

In the literature analogous transformations have been studied for the dual problem of tri- 
angulations and simplicial decomposition. In particular, it is known that, given two different 
partitions of the D-dimensional space in Z)-simplices (where a 0-simplex is a point, 1-simplex 
an edge, 2-simplex a triangle, 3-simplex a tetrahedron, etc.), one can be transformed into 
the other by a finite sequence of two local transformations called "Alexander moves" [P7|. 



The first move is the addition of a vertex inside a simplex dividing it in D -|- 1 simplices with 
the same boundary of the original simplex. The second move correspond to add a vertex 
on an edge of a simplex and connecting it with the vertices of the incident simplexes. For 
the 2D case the two Alexander move are the insertion of a new vertex inside a triangle and 
the insertion of a new vertex on an existing edge. They are shown in fig. (2a) and (2b). In 
the dual froth these moves corresponds to a cell-division which inserts a triangle near to an 
existing vertex and to a cell-division which inserts a square near to an existing edge (fig. (2c) 
and (2d)). In 3D, one can easily see, by following the same procedure illustrated for 2D, 
that the two Alexander moves can be obtained in the dual space of the froth by applying 
cell-division transformations. In the general case, one can see that the first Alexander move 
can be always done in the dual froth by dividing a cell in the proximity of a vertex inserting 
in this way a new polytope with D + 1 neighbours (a simplex). The second move can be 
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done dividing a cell in the proximity of an existing D — 1 dimensional interface between two 
cells. In this case, the kind of polytope inserted depends on the local configuration. 

We have therefore shown that, the two Alexander moves are reduced in froths to two 
special kinds of cell-division transformations. Consequently, the entire set of all the pos- 
sible froths tiling a given manifold can be generated by cell-division and its inverse (cell- 
coalescence) transformations. 

Now we investigate how the average properties of the structure are modified by these 
transformations. First consider the cell-division transformation in the two dimensional case. 
The cut of a cell corresponds to insert in the system 1 additional face, 3 edges and 2 vertices. 
Therefore, one has the transformations C2 ^ C2 + 1, Ci — Ci + 3 and Co Co + 2. One can 
verify that the Euler-Poincare characteristic rests unchanged (indeed, Xd = Co — Ci + C2). 
On the contrary, the average coordination number (which is given by X2 = see Eq.(^) 
is modified 

X^ = X2T^^iX2-6) , (8) 

with the upper sign corresponding to a cell-division transformation and the lower sign to its 
inverse (coalescence). 

Now consider the 3D case. A cell-division corresponds to insert inside a cell a new face 
which can have, in general, Ci edges and Co(= Ci) vertices. This cut corresponds to the 
transformation C3 — C3 -|- 1, C2 — >■ C2 + 1 + Ci, Ci ^ Ci -f Ci -t- Cq and Cq ^ Cq -\- Cq. 
Equation (^ gives X3 = and therefore we obtain that the average coordination number 
transforms as 

In the D-dimensional case, the cut of a cell corresponds to introduce a D — 1 dimensional 
interface which is, in general, made of Co vertices, Ci edges, C2 faces, C3 3-dimensional cells 
.... C£)_2 (-D — 2)-dimensional polytopes. Consequently, the division of a D-dimensional 
cell (or the coalescence between two cells) of the D-dimensional froth correspond to the 
transformation 

Co —>■ Cq ± Cq 

Ci Ci ± Ci ± Cq 

C2 — > C2 i C2 i Ci 

Cfc Cfc ± Cjfc ± Cfc_i 



Cd-1 Cd-1 ± 1 ± Cd~2 

Cd ^ Cd±1 , (10) 

where the upper sign (-I-) corresponds to a cell-division transformation and the lower sign 
(— ) to its inverse (coalescence). By substituting into Eq.(P one can verify that the global 
curvature {xd) is an invariant quantity under the transformation (p!0|). Note that, expression 
(PUj) takes the canonical form C^ — ^ C^ ± ± Ck-i for all k if one impose c_i = 0, c^^i = 1 
and Cd = 0. 
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From Eq.(|^) one has the identity Xj^Ck = {D + 2 — k)Ck-i (where we used the definition 
Xk = (^fc-i,fc))- By substituting in this expression the transformation (0) we get 

X'^ = XkT ^ ^ iXk{ck + Ck-i) -{D + 2- k){ck-i + Ck-2)) , (11) 

(upper sign, cell-division; lower sign, cell-coalescence). 

We recall that through cell-division/coalescence transformations it is possible to generate 
the full class of froths tiling topologically identical manifolds. The modification of the aver- 
age structural properties associated with these geometrical transformation are algebraically 
given by equation (|TT1). By using this expression it is therefore possible to find the average 
topological properties of all the froths on a given manifold. 



4 Fixed points 

When cell-division or coalescence transformations are performed on a 2D froth with average 
coordination X2 = X^ = 6, they leave the local average structural properties unchanged 
(i.e. X2 = X2 = X2 = 6, see Eq.(P)). This is a fixed point in the transformation ( |T0|) 
and corresponds to 2D Euclidean froths. Moreover, one can see that the average structural 
properties of froths which are tiling elliptic surfaces (i.e. the one with X2 < 6) are modified 
toward the Euclidean structure {X2 < X2 < Q) by the application of the cell division 
transformation. Analogously hyperbolic froths {X2 > 6) are also modified towards the 
Euclidean structure (6 < X2 < X2) (see fig. (la)). (Note that the global curvature remain 
always unchanged. Indeed, xd is an invariant under the transformation ([lOD). 

In the general case, one can immediately see that transformation (O) has the fixed point 



X; = {D + 2 -k f"-'^ , (12) 

which is the structure that is invariant under cell-division/coalescence transformations {{X^'} 

A froth is a graded set. Therefore the D — 1 dimensional interface that have been intro- 
duced into the system to cut a cell is a (D — 2)-dimensional eUiptic froth with cq vertices, 
Ci edges .... cd-2 {D — 2)-dimensional cells. All the relations written above, and in par- 
ticular Eqs.(|]) and (^), can be applied to this {D — 2)-dimensional elliptic froth. One has, 
CfeXfc = {D — k)ck-i and Ck-iXk~i = {D+l — k)ck-2, with Xk and Xk-i the average coordination 
numbers of the k and (A;— l)-dimensional polytopes which are making the (D — l)-dimensional 
interface. By substituting into Eq.(|T2D, one gets 

^^_ (D + 2-k)iD + l-k + Xk-i) 

^D + l-k){D-k + Xk) " ■ ^ ' 

The fixed point configuration is therefore determined by a set of variables {xk} with k < D 
which are the average coordinations of the {D — l)-dimensional polytope that is inserted 
or removed during the cell division or coalescence transformation. For example, in D = 3, 
relation (0) gives 

X; = 2x2 + 2 (14) 

with X2 the number of edges of the face that is inserted (or removed) to divide a cell (or 
make coalescence between two cells). 
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The minimum number of edges per cell is 3. Therefore from Eq.(|T^) follows that fixed 
point structures are possible only in the region of the parameter space with X3 > 8 (see fig.(|l] 
a)). Any structure with X3 < 8, is transformed towards the fixed point region (X3 > 8) by 
applying cell-division transformations. 

In four dimensions Eq.(|r^) gives 

^ 20X2 
' 3(2 + X2) 

12 

XI = 2(X3 + 1) = 2(- + 1). (15) 

We can express the parameter X2 in Eq.(|l5D in terms of X| obtaining XI = 5 g_^^ , which 
is the condition on the even valences that identifies the Euclidean region in 4-dimensional 
froths (see Eq.(|^)). The fixed point structures are Euclidean. Since > 4, follows X'l > 10, 
which implies that structures generated by cell division can only access to a part of the 
Euclidean region in the phase-space. 

We can in general prove that the fixed point {X^} given by Eq. ([T3|) , is the average struc- 
ture of a D-dimensional froth which is tiling a manifold with xd = 0. Indeed, let substitute 
into Eq.(|^) the fixed point configuration {{XI}) and apply the cell division transformation. 
By definition the sequence {X^} doesn't change, whereas the total number of cells increases 
of a unity {Cd ^ Cd + 1). To satisfy Eq.(|) before and after this transformation one must 
have Xd = 0. Which prove the theorem. 

By re- writing Eq.(|ll]) in the form 

Xl=X,T ^^l^"l-'^ {X,-X;) , (16) 
Cfc ± Cfc ± Ck-1 ^ ' 

it is easy to see that the fixed points are stable under cell division transformations (upper sign 



in Eq. dTq) ) which inserts identical polytopes as interface. Indeed, from Eq . (|16D , if X^ > X^ 
then Xfc > X[ > X^ and vice versa. 

In the space of the configuration { X2, X4, X21, ... }, when D even, the Euclidean region 
is a surface given by Eq.(^) (with xd = 0). The fixed point configurations are a subset of 
this surface. Froths outside the fixed point configuration are always transformed toward this 
subset by applying cell-division transformations. When D odd, Eq.(^ is not a constraint 
and the fixed point are associated with a sub- volume of the whole accessible parameter space. 



5 Construction of Euclidean froths 

In this paragraph we study froths generated by cell-division transformations. We therefore 
study the class of structures {X^} given by Eqs.(|lT]) and (|13]). The full class of these froths is 
obtained by varying in Eq.(|r^) the parameters {x^} in the allowed range (i.e. k+1 < Xk < 00, 
which satisfy the conditions (^) for k odd and the relation (^) with Xd < 0). Here, we study 
only some particular cases. 

Let us first note that, from Eq.(|T3|), the average number of neighbours of the fixed point 
structure is given in term of the coordination of the inserted interface {xd-i) by 

= 2XD-1 + 2 . (17) 
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In two dimensions an edge is inserted or removed from a face. The "coordination" of an 
edge is its number of vertices: x^^i = Xi = 2. Therefore = 6, as should be in Euchdean 
space. In three dimensions a face is inserted in, or removed from a cell. The coordination of 
this face (0:2) is its number of edges and in principle it can be any number between 3 and 00. 
But, a face with a large number of edges can be inserted only in a cell with a large number of 
neighbours and it can be removed only if it exists in the froth. Therefore, only some values 
of X2 are admissible. One can easily see that a triangle (x2 = 3) can always be inserted in 
the proximity of a vertex. Analogously, a square {x2 = 4) can also be always inserted in the 
proximity of an edge. From ([TTD follows therefore that three dimensional Euclidean froths 
with 8 < X3 < 10 can always be generated. But, in general, it should be possible to insert 
faces with higher values of X2- To have an estimation for the "typical" value for the number 
of edges of the inserted face let make a cut of the whole three dimensional froth with a plane. 
The result is a two dimensional Euclidean froth where each single face is the result of a cut on 
a three dimensional cell. This two dimensional froth is therefore a representative set of faces 
produced by random cuts of three dimensional cells. The average number of edges for this 
set of faces is X2 = = 6. Therefore, from (0), a "typical" froth generated by cell division 
is expected to have a fixed point coordination around = 14 p8[. Cells in biological 



tissues appear in various polyhedral shapes with a number of faces distributed in a narrow 



range around 14 |29|. A widely studied three dimensional froth made with identical cells 



is the "Kelvin froth", its cells are space-filling truncated octahedra with X3 = 14 pO| , pT 
Coordinations between 15.53 and 14 are found in Voronoi partition of space where the 
higher value corresponds to a Voronoi partition from random points |3^ whereas the lower 



value corresponds to more compact and homogeneous packings . Smaller values in the range 
13.333 < X3 < 13.5 characterize an interesting class of natural structures (Frank-Kasper 
phases [^) which partitionate the ordinary space with cells with pentagonal and hexagonal 
faces only. Soap froth has typically X3 ~ 13.7 P^ . 



In a D-dimensional froth a — l)-dimensional interface with coordination x^^i is in- 
serted or removed by cell division or coalescence transformations. As pointed out above, 
simplexes with coordinations xd-i = D can always be inserted in the proximity of an exist- 
ing vertex. This is the minimum possible value for xd-i and substituted into Eq. ([T7|) sets 
the minimum value of the average number of neighbours in a D-dimensional Euclidean fixed 
point structure at the minimal coordination 

X}) = 2D + 2 . (18) 

The argument for the "typical" cut that we used in three dimensions can be directly 
extended to any dimension. Indeed, one of the properties of froths is that a cut with an 
hyper-plane of a D-dimensional froth generates a. D — 1-dimensional Euclidean froth. For 
instance, a cut of a four dimensional froth give a three dimensional Euclidean froth. We can 
assume, that this froth has the "typical" coordination X^ = 14 found above. Inserting into 
Eq.(|17D one gets X^ = 30. The same arguments, extended to any dimension, give 

X*^ = 2^^+^) - 2 , (19) 

for the average number of neighbours per cell in the "typical" D-dimensional froth. 

What makes Eq.([T3|) powerful is the fact that, not only the average number of neighbours, 
but all the average properties of the fixed point structures can be deduced in term of the 
properties of the inserted interface. 
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5.1 Minimally coordinated froths 

Let first construct the Euclidean froth with minimum coordination numbers. It is the fixed 
point structure associated to a cell-division transformation which inserts interfaces with 
minimal coordinations. This interfaces are D — 1 dimensional simplices inserted in the 
proximity of vertex. They have Xk = k + 1 (for k < D). By substituting into Eq.(0) one 
obtains 

(Note that = 2D + 2, as discussed above.) This are the average structural properties 
of a froth which is tiling a manifold with xd = which is homologue to the Euclidean 
space. It is the known Euclidean forth with minimal coordination numbers. Starting from 
any given D-dimensional froth one can always transform it into this minimally coordinated 
one by applying an infinite number of cell-divisions near existing vertices. The resulting 
structure is expected to have cells with very different topological properties. Indeed, for each 
cell-division transformation a new cell with D + 1 neighbours is inserted and 1 neighbour is 
added to the D + 1 cells around the inserted simplex, distributing therefore the coordinations 
inhomogeneously between cells. 

5.2 Homogeneous partitions 

A D-dimensional froth has -D + 1 edges incident on each vertex. In an ideal homogeneous 
partition of space these edges are equally separated in angle. That corresponds to an angle 
Qideai _ cos~^(— 1/D) between each couple of edges |2^. In a froth, edges must close in rings 



which are bounding two dimensional faces. It is easy to see that with the angle Q'^'^'^"-^^ fiat 
rings close with an average number of edges equal to: 

X't'' = ■ (21) 

vr - COS n - 

This number is 6 in two dimension, 5.104 in three dimensions, 4.767 for D = 4, 4.588 for 
D = 5 and tends to 4 when D ^ (yo. Note that Xlf^"-^ is irrational for D > 2. In the 
Euclidean space, the "ideal" structure cannot be obtained by any ordered lattice structure. 
This is an example of geometrical frustration. But disordered or non-periodic structures can 
approximate X2^'^°''' with arbitrary precision avoiding in this way the frustration. 

The average number of edges per cell Xlf^"-^ (Eq.(|21|)) is the only quantity, of the ideal 
structure, that can be calculated by using these geometrical arguments. All the other coor- 
dinations are unknown, but we can construct fixed point structures that approximate this 
ideal froth in the ring coordination X2. For these structures we can calculate the whole set 
of coordinations, and therefore we can infer information about the coordinations of the ideal 
one. 

We expect that structures that partitionate uniformly space must have X2 — X^^°'K 
By cell division transformation it is possible to generate froths that approximate the ideal 
structures by inserting interfaces with X2 close to the ideal value Xlf'^'^K 

For = 3, Xlf^"-^ = 5.104 which corresponds to Xg*^"^"' = 13.392. In 3D we can generate 
homogeneous partitions by inserting pentagons 0:2 = 5 or hexagons X2 = Q obtaining (from 
Eq.(|l7D) 12 < X3 < 14, which is in the right range. 
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In general, since we are looking for homogeneity, it is logical to insert as interfaces regular 
polytopes with X2 close to the ideal value Xlf^"-K These polytopes can only be hyper-cubes 
{4,3...} (which have X2 = 4) and the polytope {5,3...} (with X2 = 5), but this second polytope 
exists only up to D = 4 ([^ and App. |^). 

(a) Cell-division operations which inserts the polytope {5,3...} can therefore generate Eu- 
clidean fixed point structures up to D = 5. In D = 3 this corresponds to cell-divisions which 
insert pentagonal faces, that generates a fixed point structure with = 5 and X^ = 12. In 
four dimensions, the fixed point structure obtained by inserting dodecahedra ({5,3}, X2 = 5, 
X3 = 12) has X2 = 100/21 = 4.7619... value that is very close to the ideal one. The four 
dimensional cells of this froth have X^ = 26 neighbours in average. For D = 5, hj dividing 
cells with the polytope {5,3,3 } {x2 = 5, X3 = 12, X4 = 120) we obtain a fixed point structure 
with X2 = 75/16 = 4.68... (see Eq.(|T3|)) which is larger than the ideal value. The average 
number of neighbours is in this case X^ = 242. 

(b) The average coordinations of the fixed point froths generated by inserting hyper-cubes 
{4,3,...} are given by imposing Xk = 2k into Eq.(|2^)) 

^ ' ^ {D + l-k){D + k) ^ ' 

In this structure the D-dimensional cells have = AD — 2 neighbours in average. The 
average ring-coordination is = 4:D{D + 1)[(-D — 1){D + 2)]^^ which correctly tends to 4 
when D ^ 00, but it is systematically lower than X2^'^°''' for D > 2. This is presumably a 
rather inhomogeneous structure 

(c) To maximize homogeneity one can construct a structure by inserting new interfaces with 
the same topological properties of the existing structure. We expect a resulting structure 
that evolves toward a self-uniform homogeneous partition. Let therefore perform cell-division 
transformations by inserting interfaces with Xk = X^ (with k = 2, D — 1). By substituting 
into Eq. (|^) one obtains a recursive equation with the following solution 

(c) ^: = (l + 7J^^)^ • (23) 

Here Xa = D(D + 1) and X2 = 2^^^, which asymptotically tends to 4 and is much closer 
to the ideal value than the one of structure (b). 

(d) Partitions can be generated by inserting "typical" interfaces as described before. Here 
the "typical" {D — l)-dimensional interface has coordinations Xk which are equal to average 
coordinations of the fixed point Euclidean structure obtained with this procedure in D — 1 
dimensions. The values X^ are then given in term of a recursive equation (with initial 
condition X^ = 2). Here are the solutions for k = 2 and k = D, which have a simple 
compact form 

(d) x; = 12 ^ 



3D -2 

XI = 2(^+^)-2 . (24) 
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Surprisingly the product of these valences from k = 1 to D has also a very simple form: 
Xj'XIXI ■ ■ -X^ = D\{D + Here, the value of X2 is systematically bigger than the one 
of the ideal structure but it is extremely close to it. 

In table 1 the values of X2 and X^ are reported, up to = 10, for the whole set of 
fixed point structures which have been studied in this paragraph. In Fig. (^) the value of X2 
for the ideal partition and for the fixed point structures (a), (b), (c), (d) are plotted up to 
D = 24. 



5.3 Kissing Numbers 

The average number of neighbours X^ of a D-dimensional cell in an Euclidean froth is an 
interesting quantity. In sphere packings a corresponding quantity is called "kissing number" , 
it is the number of identical spheres that can be placed around a given sphere being in contact 
(kissing) with it [jl3|. To find sphere packings configurations with high kissing numbers has 



relevance in the design of efficient codes. It is known that, for packings of identical spheres, 
the kissing number (KN) is 6 in D = 2, 12 in D = 3, but exact answers are unknown 
for dimensions above 3 except for D = 8 (KN=240) and D = 24 (KN=196560) where two 
specially dense lattices Eg and achieve the maximal possible values of KN. In fig.(|) 
are reported the values of the highest known kissing numbers for lattice and non-lattice 
sphere packings. Two known bounds for KN when D ^ 1 are also reported. The kissing 
number question concerns to find the best local arrangements of spheres. In high dimensional 
spaces, this configuration does not necessarily corresponds to any lattice packing. Disordered 
or quasi-ordered packings are often more suitable to attain high kissing numbers. Dimension 
D = 9 is the first where non-lattice packings are known to be superior. Here the Leech 



lattice Ag has KN=272 whereas the best bound known is 380 ||13 . 

To any sphere packing one can associate a cellular structure constructed by partitioning 
the space in convex polytopes each one containing inside a sphere. Kissing spheres are 
neighbours. In a dense sphere packing the enveloping polytopes make a space-filling partition 
of space. The number of neighbours of this system of polytopes is related with the kissing 
number and it is expected to be bigger than KN, because some non-kissing spheres can 
be first neighbours in the associated froth. This is for instance the case in D = 3 where 
the configuration with KN=12 corresponds to a close packing of spheres with an associated 
Wigner Seitz cell that do not pack in a froth: the incidence numbers are not minimal. 
This is a topologically unstable configuration. Infinitesimal random displacements change 
the number of topological neighbours from 12 to an average value of 14, but in this case 
neighbouring spheres will be not all in contact. In general, in close packings, we expect the 
number of neighbours X^ of the enveloping polytopes to be bigger, but of the same order of 
magnitude, of the kissing numbers of the enveloped spheres. 

In fig.(^ the kissing numbers for some known sphere packings are compared with the 
coordination numbers obtained from our homogeneous partitions (a), (b), (c) and (d) up to 
D = 24. 
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5.4 Voronoi partitions 



The average number of vertices on the boundary of a D-dimensional cell {(rio^o)) can be 
exactly calculated for Voronoi partitions generated from random points ^ : 



(^o,d) 



r{D) 



D 



T UD + r 



r i r iD + i 



r(i(z^ + i; 



D-l 



(25) 



Asymptotically this quantity scales as (?t.o,d) oc D^^'^ ^ 

The average number of vertices per cell can be expressed in term of the coordinations by 



(26) 



(Note that for D = 3 Eq.(^) give (t^o.d) = 967r^/35 that substituted into Eq.( ^6]) leads to 
X3 = 2 + 487r^/35 = 15.53...). By substituting the fixed point configuration into (^61) we find 
that the structure (b) has {uo^d) = {D + l)2^~^, the structure (c) gives (n^D) = (2-D)!(-D!)^^, 
whereas (d) has (/io,z)) = (-D + 1)!. In figure (§) the behaviors of the average number of 
vertices per cell in the Voronoi froth and in the three structures (b), (c) and (d) is shown 
for 3 < D < 50. 



6 Conclusions 

The topological structure of a D-dimensional cellular system can be characterized, in average, 
in terms of the coordinations (Xk) of the irregular polytopes which are making the structure 
(Section Only the coordinations of the even dimensional polytopes (^2^) are necessary 
for this characterization, the odd ones are expressed in terms of the even ones by the relation 
(^). Therefore, the average structure of a D-dimensional froth is characterized by a sequence 

{X2, Xi, } of Y — 1 (or -^^) variables for D even (or odd). These variables are 

related with the space-curvature through Eq.(^. Regions in the parameter space {X21} 
corresponding to D-dimensional froths tiling spaces of different curvature are discussed for 
D < 4 (figdU) and Eqs.(D and (0)). 

We use cell division and coalescence transformations to build D-dimensional froths (Sec- 
tion]^). We show that through cell-division/coalescence transformations it is possible to gen- 
erate the entire class of froths tiling topologically identical manifolds. The dynamical renor- 
malization of the variables X^. under such transformations is found (Eq.([lID). The existence 
of classes of structures which are invariant (fixed points) under cell division/coalescence is 
pointed out (Eqs.(|T^ and (|13|) ). We show that these structures are tiling Euclidean spaces. 

Several fixed point Euclidean structures are constructed in Section ^ We discuss the 
average statistical properties for the of minimally coordinated Euclidean froths, and for 
several topologically homogeneous space partitions (Eqs( PIl| - P^ ) and tab.([^)). The topological 
properties of the most homogeneous cellular partition are searched and compared with known 
geometrical results (fig.(|^)). 

Finally, the fixed point Euclidean structures are compared with known high dimensional 
structures generated by sphere packings and Voronoi constructions (figs.(|) and (|^)). 
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A The existence of polytopes {5,3, ...} up to D = 4 

There is a class of regular polytopes with X2 = 5 and minimally connected vertices ({5, 3, ...}) 



which exists up to dimension D = 4 |26 



They are pentagons in D = 2, dodecahedra ({5, 3}) in D = 3 and polytopes {5, 3, 3} 
in D = 4. A tessellation of pentagons makes a 2D elliptic froth with X2 = 5 and C2 = 
12, it is a dodecahedron {5,3}. A tessellation of dodecahedra make a 3D elliptic froths 
with X2 = 5, X3 = 12, C3 = 120, which is the {5, 3, 3} structure. It turns out that a 
tessellation with {5,3,3} polytopes does not makes any 5-dimensional polytope [^. If 
existing, such a structure would be a 4-dimensional polytope with X2 = 5, A3 = 12, A4 = 
120. By substituting these values into Eq.(|^), we get C4 = X4; which implies Xi > 0. This 
hypothetical structure would therefore be an elliptic froth, homotrope to a sphere which 
implies Xi = 2. Then, from the previous identity, C4 = 2. The hypothetical {5, 3, 3, 3} 
structure would be an elliptic froths which closes onto itself with two cells only. But two 
cells are insufficient to make a 5-dimensional polytope (the minimum number is 6). It follows 
therefore that the {5, 3, 3, 3} structure does not exists and neither exist the others higher 
dimensional tessellations {5, 3, 3, 3...}. 
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Table 1: Average ring coordination {X2) and cell coordination {Xd) for some Euclidean fixed 
point structures (FP) generated by cell division (see text). 
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Exislence of 3D froLli 




Figure 1: Different regions in the parameter space {X21} corresponds to froths which are tihng spaces 
of different curvatures, (a) two dimensional froths with X2 < 6 tile elliptic surfaces, froths with X2 > 6 
tile hyperbolic surfaces and X2 = 6 corresponds to froths tiling the Euclidean plane, (b) The hyperbolic, 
Euclidean and elliptic tilings corresponds to three regions of the {X2,X4} parameter space. Cell division 
transformations modify the properties of curved tilings towards the Euclidean ones (arrows). 
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Alexander moves 




Complex Froth 



Figure 2: Given two different partitions of the _D-dimensional space in D-simplices one can be transformed 
into the other by a finite sequence of two local transformations called "Alexander moves" ((a) and (b) for 
the two dimensional case). In the dual froths these transformations corresponds to two special cell division 
transformations ((c) and (d) for the two dimensional case) 
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Figure 3: The average ring connectivity X2. for homogeneous partitions of spaces of dimension D, obtained 
from a geometrical approach [25] (full line) is compared with the one of the fixed point structures (symbols). 
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Figure 4: Kissing Numbers KN (maximum mimber of identical hard spheres that can stay around and touch 
a given sphere in D dimensions) and coordination numbers Xd (average number of cells around a given cell 
in a D dimensional froth) are compared for some known sphere packings (open symbols) and for the fixed 
point structures (full symbols). Two upper and lower bounds for the KN in high dimensions are also plotted 
(full and dotted lines). 
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Figure 5: The average number of vertices in a D dimensional Voronoi' cell generated from Poissonian points 
can be exactly calculated [32] (black circles). Here this number is compared with the one associated with 
the fixed point structures (gray symbols). 
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